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We identify the Goldstone modes appropriate to the low and high density phases 

of SU{2) lattice gauge theory with staggered fermions. We present hybrid Monte 

f^ ' Carlo simulation results for susceptibilities on a 6* lattice at /3 = 1.5, m = 0.05. 

Q^ ' We specify how to implement a diquark source term in a lattice simulation and 

[j^ , present first measurements of the lattice diquark condensate as a function of the 

, diquark source. 



The symmetry breaking pattern of SU(2) lattice gauge theory with fermions 
in the fundamental representation of the group is different from that of the con- 
tinuum model where the symmetry breaking nattern for fermions in a pseu- 
(^ ' doreal representation is SU{2N-f) -^ Sp{2Nf)u. To identify the symmetries 

Cn , appropriate to the lattice model with staggered fermions we start with the 

^^ ' kinetic term in the lattice action for a gauged isospinor doublet □ of staggered 

^i , fermions. For clarity we consider iV = 1 flavors. 
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Skin= ^ -^^^ [x{x)U^{x)xix + i>) - x{x)Ul{x ~ i>)x{x ~ £>)] 

X. z-*— 1,3 



+ E^ x{^)e'Utix)x{x + i)-xix)e-^'U}{x-i)x{x-i) .(1) 



^ ■ Defining: 

C^: X.^ixe,-xTr2) ■■ X,= ( ^° ) (2) 

and using r]^{x ± fi) = rj^{x) and T2U^T2 = U* (where T2 is a Pauli matrix) 
the kinetic term can be expressed in the following form: 
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Xe{x)(^^^ f_^jUt{x)Xo{x + i)- 
Xe{x)(^^l' g° ) Uj [X - i)Xo {X - i) 
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There are two manifest global U{1) symmetries in the original action (pi): 

X I— > e*"x X '"*■ X6~'" '■ X '"^ e*"'^x x '"^ X^*"^ (4) 

where e{x) — (—1)^1+^2+2:3+3:4 rpj^g ^^.g^ invariance corresponds to baryon 
number conservation: the second which holds only in the chiral limit (m — > 0) 
corresponds to conservation of axial charge. Both of these are subsumed in a 
larger U{2) symmetry which only holds for m — fj, = 0: 

Xo^VXo Xe^XeV^ V"eU(2). (5) 

In a continuum approach the axial anomaly would reduce the symmetry to 
SUi2). 

For ^ 7^ the full U{2) lattice symmetry is reduced to U{1)v(^{^)a and is 
further reduced to U{l)v for /i 7^ 0, m 7^ 0. It has been proposed^ that U{l)v 
is likely to be broken spontaneously by the formation of a diquark condensate. 
At fi = the model displays spontaneous chiral symmetry breaking, with a 
chiral condensate of same form as the mass term ^^ x{^)x{x)- 
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Xe( ^ ^ JT2XT+XT[^ ^ \T2X0 
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This chiral condensate breaks the global symmetry (^ . The residual sym- 
metry is generated by the subgroup which leaves J invariant of which 



the most general element is I _jq, j, which generates a C/(l). Thus we 

identify the pattern of chiral symmetry breaking appropriate to the lattice 
model with iV = 1 as f7(2) — > V(X). Now our hybrid Monte Carlo algorithm 
simulates with det MUf and since detAft = detA/* = detT2AfT2 = detAf, 
this describes two identical staggered fermion flavours resulting in pattern of 
symmetry breaking C/(4) -^ 0(4), with 10 associated Goldstone modes. 

At high density, we postulate that a large Fermi surface will promote the 
formation of a diquark condensate. In principle many diquark wavefunctions 
can be written down therefore it is a dynamical question as to which conden- 
sate actually forms. We can deduce the "maximally attractive channel" for 
the diquark condensate by insistincuthat the condensate wavefunction is anti- 
symmetric under exchange of fieldsQ and by assuming the condensate is gauge 
invariant, invariant under lattice parity and as local as possible in the x fields. 
A local condensate qq-2 which satisfies these conditions is: 

mi = 2 U*"(a^)^2x(a;) + x(a;)T2X*"(a;)] . (7) 



In terms of the X, X fields the diquark condensate is written 
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Xe[ ^ At2XT+XI^(^ _^ ]t,Xo 



(8) 



The diquark condensate (Jg) can be obtained from the chiral condensate 
(^ by an explicit global U(2) rotation: 

Therefore in the limit m = 0, /x = the pattern of symmetry breaking remains 
U{2) -^ U{1): as /i increases the condensate simply rotates from a chiral qq to 
a diquark qq. 

The full set of Goldstone modes (again for the simplest case TV = 1) can be 
found by considering infinitesimal rotations of either the chiral condensate (ph 
in the low density phase or the diquark condensate (0) in the high density phase 
by V^ = 1 + i6X, with A one of the U{2) generators {l,Ti}, and identifying 
the mode as the coefficient of 0{6). In the chiral limit and for /i — we 
expect 3 massless Goldstone modes (since dim U(2)=4). Since the full U{2) 
symmetry of the action is reduced for ^ ^ only two of the four generators 
correspond to Goldstone modes i.e. the generators appropriate to U{l)v and 
U{1)a which are T3 and 1 respectively. In the low density phase the rotation 
generated by 1 gives the pion: x^X: while the rotation generated by T3 leaves 
the condensate invariant. Rotation of the diquark condensate by 1 gives a 
pseudoscalar diquark: x^^T2£X + X'''2£X*'^jWhich is pseudo-Goldstone for fi ^ 
and rotation by T3 gives a scalar diquark: x''^T2X — X'''2X*'^ ? which remains an 
exact Goldstone mode even for to 7^ 0. The fact that the low and high density 
regimes are characterised by different numbers of massless modes indicates 
that they should be separated by a true phase transition at some 11 — fic- In 
general we expect fic = "mib/Nc where mib is the mass of the lightest baryon 
and Nc is the number of colours. The SU{2) theory is a special case because 
we expect fic — titt/S i.e. there is a Goldstone baryon in the theory due to 
the SU(2) symmetry which relates quark to anti-quark. In this respect the 
model is crucially different from QCD where we expect fic — Triproton/'^ yet 
find fic ~ ^7i7r/2 due to pathologies in the hybrid Monte Carlo simulations for 

We can formally introduce a diquark source term into the fermionic sector 
of the partition function as follows: 

^form. = / dipdip exp {i;Mip + ipJtp + -ipJip) (10) 
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where M is the fermion matrix while J and J are diquark and anti-diquark 
source terms respectively. It is necessary to recast this equation in order to 
perform the integration over the fermion fields. To do so we define a new 2 



component Grassmann field: 

terms of a "double matrix" 



and then reconstruct Zform. in 



Zfe 



d(f){4' ip) 



J M. 

2 



M 
2 

J 



dot 



jr. f 
-^ jr. 



1/2 



(11) 



where we have used J d(pexp((pA(j>) ~ vdctA = Pfaflian(A). Note that we 
require an antisymmetric matrix therefore we chose J = J = JT2 with j a real 
number. 

Making use of the property Indet A — Tr In A, on the double matrix leads 
to the lattice diquark condensate: 
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In Fig. 1 we show the diquark condensate as a function oi j and /i. Polyno- 
mial fits to these curves were used to obtain the j ^ limit plotted alongside 
the standard observables in Fig. 2. The sharp rise in the diquark condensate 
at /i ~ 0.4 coincides with the downward jump in the average plaquette while 
the chiral condensate decreases smoothly as /i is increased. 

Measurements of the susceptibilities (Figs. 3,4) are strongly suggestive 
of a phase transition at fic = 0.4. This is consistent with the conclusion 
from the observables in Fig. 2. We shall refer to /i < ?Ti7r/2 ~ 0.4 as the 
low density phase and /x > TO,r/2 as the high density phase. We find, as 
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expected, that at /i = the scalar diquark is degenerate with the pion while the 
pseudoscalar diquark is degenerate with the scalar meson. The pscudoscalar 
diquark which is heavy in the low density phase becomes light in the high 
density phase whereas the pion which is a pseudo-Goldstonc in the low density 
phase becomes heavy in dense (chirally symmetric) phase. As seen in Fig. 
4 susceptibilities of scalar and pseudoscalar diquarks become very large in 
the dense phase. In the dense phase we observe that the scalar diquark and 
the pseudoscalar diquark are light states while the connected contributions to 
the scalar imply that it is comparatively heavy. In our symmetries analysis 
we identified the scalar diquark as an exactly massless Goldstone mode in 
the dense phase while we expected the pscudoscalar diquark to be a pseudo- 
Goldstone mode for m 7^ 0. The scalar meson is not expected to be a Goldstone 
mode in the high density phase. Although the pion and scalar meson appear 
to be approximately degenerate in the high density phase this degeneracy may 
be broken when the disconnected contributions to the scalar are taken into 
account. 
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